"de Gennes" narrowing in supercooled molecular liquids : Evidence for 
center— of— mass dominated slow dynamics 
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The density correlation function <J>(q, t) of the two similar substituted aromatic liquids, Toluene 
and m-Toluidine, is studied by coherent neutron spin-echo and time-of-flight scattering for wave 
vectors q around the maximum g max of the total static structure factor Sm(q) in the supercooled 
region. The wave-vector dependence of the mean structural a-relaxation time (r 9 ) shows in both 
liquids a very pronounced de Gennes-like narrowing centered around go < g ma x, where go corre- 
sponds to the first maximum in the center-of-mass static structure factor Scom(o). We find that 
the narrowing can be described quantitatively by using ScoM(q)/q 2 instead of 5' m (g)/g 2 indicating 
that at the corresponding molecular length scales the relaxation of 4>(g, t) is dominated by purely 
translational motion. 
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It is well known that the half width of the experimen- 
tally observed dynamic structure factor S(q, u>) (wave 
vector q and frequency to) of classical dense monoatomic 
fluids in equilibrium shows a pronounced minimum at 
wave vectors q ~ g max w 2ir/a (a being the atomic di- 
ameter) where the static structure factor S(q) has its first 
maximum. This effect, called "de Gennes narrowing" [Q, 
was suggested by de Gennes in 1959 on the basis of sum- 
rule arguments. It implies that plane- wave density fluc- 
tuations with q in the vicinity of g max decay slower in 
time due to the strong spatial correlations that exist at 
these wave vectors. Experimentally, this effect has been 
observed by coherent quasielastic neutron scattering on 
simple atomic fluids like noble gases Q and also on liquid 
metals || . It has been confirmed by computer simulations 
on simple fluid models 

The extension of this simple concept to molecular liq- 
uids in the supercooled, i. e. high-density state, is difficult 
and far less investigated . Experimentally, only a few 
studies of the wave vector dependence of the a-relaxation 
in supercooled liquids indicate that a narrowing is 
present. The effect, however, was small and the analy- 
sis remained only qualitative. Moreover, in glass forming 
polymers no evidence for any clear modulation was re- 
ported in an extended analysis where all relaxation pro- 
cesses were taken into account . The difficulties mainly 
arise from the following reasons: (i) in the supercooled 
regime, the presence of memory effects leads to a struc- 
tural a-relaxation with a marked stretching of the corre- 
lation function; (ii) for molecular liquids the intermolecu- 
lar structure factor S m (q) is a linear combination of many 
partial atom-atom structure factors; (iii) the dynamical 
observables depend in general upon both translational 
and rotational degrees of freedom. 

In this letter we present the first quantitative coher- 



ent neutron scattering study of two very similar sub- 
stituted aromatic liquids Toluene and m-Toluidine. A 
large g-range was covered including the main peak of 
the static structure factor S m (q) whose maximum is at 
<Zmax = 1.8 A" 1 . 

In both liquids we find an extremely pronounced "de 
Gennes narrowing" , not yet reported for the a-relaxation 
of any glass forming liquid. It is centered around a wave 
vector qo < g max , where qo is equal to the first maximum 
of the center-of-mass static structure factor ScoM(q)- 
Moreover, the wave- vector dependence of the mean struc- 
tural relaxation time (r q ) can be quantitatively described 
by (r q ) oc 5com(?) / q 2 ■ This indicates that around these 
wave-vectors the translational motion dominates the re- 
laxation of t). 

Fully deuterated Toluene (C 7 D 8 , T 9 =117K, 
T.,„=178 K) was obtained from Euriso-top and sealed into 
quartz capillaries (Fa. Hilgenberg, Malzfeld). Optical in- 
spection showed no signs of crystallization when cooling 
down to T g . Fully deuterated m-Toluidine (C7D 8 ND 2 , 
T 9 =183.5K, T m = 243.9 K) was filled into a thin hollow 
aluminium cylinder. Note that for temperatures around 
the calorimetric glass transition the time scale of the 
methyl group rotation around their C3-axis is several or- 
ders of magnitude faster than the relevant a-relaxation 
||. Moreover, the quasielastic contribution of the rota- 
tional motion is negligible ||. Both liquids are simple 
molecular glass formers and particularly interesting as 
they can be thus considered as rigid planar and small 
organic molecules. Their dynamic has already been stud- 
ied by a large variety of experimental methods, e. g. 
viscosimetry, dielectric relaxation, nuclear magnetic res- 
onance, light- and incoherent neutron scattering ]Io|| . 

The main experiments were performed on the neutron 
spin echo spectrometer (NSE) IN11 of the Institut Laue 
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Langevin in Grenoble, France. IN11 was operated in a 
double spin echo configuration with the new multidetec- 
tor covering 30° . Different mean incident wavelengths Xi 
and scattering angles 29 were chosen to cover the q region 
around the first diffraction peak The polariza- 

tion was measured for several temperatures. Normaliza- 
tion was achieved by division by the resolution function 
measured on a quartz rod of similar shape and by the po- 
larization at t = 0. We thus obtain the normalized density 
correlation function &(q,t) = S(q, t)/S(q, 0). Additional 
experiments have been performed for Toluene on the 
time-of-flight spectrometer (TOF) Mibemol of the Lab- 
oratoire Leon Brillouin in Saclay, France. The dynamic 
structure factor S(q, u>) was measured with an elastic en- 
ergy resolution AE=36 /xeV giving access to a q range for 
elastic scattering of 0.2 < q e \ < 1.7 A -1 . After standard 
corrections the data were interpolated to constant q with 
step Aq = 0.05 A -1 to yield the dynamic structure factor 
S(q,u). 

Spin-echo data are interpreted by using the stretched 
exponential function of Kohlrausch 

$(g,t;T) = /,exp(-(t/T,) /J «), < q < 1 (1) 

where in principle all parameters can be temperature and 
q dependent. As (r g ) and (3 q are strongly correlated, we 
use the mean a-relaxation time (r q ) = / dt&(q, t)/f q — 
(3~ l T{(3~ l ){T q ) as a fitting parameter instead. 

The individual correlation functions <$>(q,t;T) shown 
in Fig. 1 indicate already without any fitting that the re- 
laxation is slowed down by about a factor of ten around 
q ~ 1.3 A" 1 when compared to other q. Due to the 
stretched character and the strong ^-dependence of (r q ) 
a much larger dynamic time window would be needed 
in order to cover all g-dependent t;T) at one given 
temperature. Although a time range of about 2.5 decades 
is covered our data do not really allow for reliable three 
parameter fits to the individual <&(q, t; T). 

However, quantitative information and consistent re- 
sults with acceptable error bars can be extracted by com- 
bining measurements at different temperatures. With- 
out altering the interpretation of the underlying physics 
we apply the well-known time-temperature superposi- 
tion principle which means that the correlation functions 
<&(q,t;T) can be superimposed simply by rescaling the 
time axis. Its operational validity at least in a narrow 
temperature range of about 30 K in the time range of 
neutron scattering was demonstrated for a large number 
of glass forming substances We use a compila- 

tion of a-relaxation times r a obtained by other spectro- 
scopic techniques [jl(| to rescale the experimental times 
according to t — ir a (To) / T a (T) , with an arbitrary nor- 
malization at T = 140 K and T = 230 K for Toluene 
and m-Toluidine, respectively. Indeed, for both liquids all 
<f>(q,t;T) collaps onto a temperature-independent mas- 
ter curve $(q,t) without any adjustment of f q (T) as il- 
lustrated in Fig. 2. The successful scaling implies that the 



microscopic relaxation time measured by neutron scatter- 
ing is proportional to the " macroscopic" relaxation time 
and that the exponents f3 q and the amplitudes f q cannot 
depend strongly on temperature. 

^From a fit with ([!]) to the master curves $>(q, t) we ob- 
tain the relaxation time (r 9 (To)) which are displayed in 
Fig. 3 normalized to its maximum value at qo- For both 
liquids, the relaxation time (r g ) has an overall variation 
with q (far from the hydrodynamic regime on the low- 
q side and from the free-particle behavior on the bigh-g 
side) of a factor of ten. It passes through a maximum as 
q is varied from 0.8 to 2.0 A" 1 . To our knowledge this 
variation is much stronger than that in any other liquid 
studied so far. More strikingly, the strength of the varia- 
tion of (t,j) with q and the position of its maximum can 
quantitatively be described by (r q ) oc <Scom(<7)/<7 2 where 
the center-of-mass static structure factor Scom(<z) was 
obtained by Monte-Carlo simulations jll],[l2|. It was not 
expected that de Gennes's very general concept jlj would 
strickly apply to molecular liquids like Toluene or m- 
Toluidine. It is far from being trivial because the corre- 
lation functions are stretched and the relaxation times 
are long compared to those of atomic fluids. In fact, 
our findings indicate that in the q-range corresponding 
to molecular length scales (van der Waals radius) and 
in the a-relaxation regime, neutron scattering mainly 
probes the motion of the centers of mass, i.e. mainly 
probes the translational motion of the molecules. This 
result is of fundamental importance for the interpreta- 
tion of neutron scattering spectra in general. What we 
actually measure are, of course, atom-atom correlations 
involving the H,D and C nuclei. In addition to a transla- 
tional component, these correlations also contain contri- 
butions involving the rotational motions. Our joint use 
of Monte-Carlo simulations and neutron scattering ex- 
periments indicates that at the q values considered and 
for the a-relaxation, these latter contributions are small- 
compared to the translational one. This finding is also of 
particular interest for microscopic glass transition the- 
ories like the mode-coupling (MC) theory, either in its 
idealized form flB]] or in the molecular version in- 
deed in MC theory the static structure factor is the only 
major input that determines the slowing down of the a— 
relaxationThe MC theory gives at least a qualitative de- 
scription of the q-dependence of (r 9 ) around the static 
structure factor maximum | fL7| and the importance of 
the center-of-mass slow dynamics in this framework is 
illustrated by the so-called semi schematic model |l8j . 

To underline the robustness of our NSE results (r q ) has 
been independently obtained for Toluene from Mibemol 
experiments by fitting ([!]) the S(q,u>) at 200 K in fre- 
quency space convoluted with the measured resolution 
function (see Fig. 3). The stretching (3 q was kept fixed 
to 0.5 close to the average value 0.48 obtained from our 
spin-echo experiment (see below). As can be judged from 
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Fig. 3, excellent agreement within a few percent between 
NSE and TOF is found for the g-dependence of the re- 
laxation time over a large q range. To emphasis this good 
agreement even further we included in Fig. 2 <fr(q,t) ob- 
tained from Mibemol spectra S(q,u>) by Fourier decon- 
volution with the measured resolution function. Thus, an 
identical wave-vector dependence of (r q ) is observed from 
the stable liquid state, where the typical relaxation times 
are of the order of picoseconds, to the supercooled state, 
with relaxation times of the order of nanoseconds. Such a 
behavior was also observed in orthoterphenyl (OTP) 

Another important output of our study is the (in- 
dependence of the amplitude f q of the a-relaxation. In 
Fig. 4, f q shows a maximum at q — 1.4 A -1 which is 
very close to the maximum of the center-of-mass static 
structure factor 5'com('?), in agreement with the numer- 
ical solutions of mode-coupling equations Ipjf . However, 
contrary to this latter, it does not oscillate in phase with 
ScoM(q) (nor S m (q)), and only a very smoothly decreas- 
ing behavior is observed at higher q due to additional 
rotational contributions |isj] . We could not detect any 
systematic g-dependence of the stretching parameter f3 q 
which is at variance to theoretical predictions |l7j and 
to the results in OTP ||. Note however, that in OTP 
no variation of (3 q with q was detected at wave num- 
bers corresponding to the center-of-mass correlations 
but at larger q values. ^From NSE we find an average 
(f3 q ) q ~ 0.48 ± 0.02 for Toluene and {(3 q ) q ~ 0.58 ± 0.02 
for m-Toluidine. The conjecture of q dependent oscilla- 
tions of P q and (r 9 ) indicates a systematic error source 
in NSE experiments: since we have used a wave vector 
band Aq/q — AA/A ~ 10% which is nearly as broad as 
the diffraction peak itself 0,^2| we have actually mea- 
sured averages (<&(q,t;T)) q of correlation functions with 
possibly different parameters P q and {r q }. In practice, a 
sum of different stretched exponentials can again be fit- 
ted by a Kohlrausch with a smaller effective exponent /3 c ff 
and a different (r 9 } c ff. ^From numerical examples using 
(t 9 ) cx q~ 2 ScoM and 0.4 < (3 q < 0.5 as input we conclude 
that the difference between (T q ) e s and (r 9 ) is always less 
than 10 % and the difference in the line-shape is anything 
between 0.01 and 0.05 which is far less than our experi- 
mental accuracy i. e. our conclusion remain valid. 

The main result of our communication concerns the 
wave vector dependence of the a relaxation time (r g ). 
The maximum in (r 9 ) at qo < g max can be interpreted 
as a de Gennes narrowing for molecular glass forming 
liquids resulting from the strong spatial correlations be- 
tween the centers-of-mass of the molecules. The relax- 
ation time characteristic is quantitatively described by 
r q oc Scom{q)/<1 2 ■ Scom(i) is the spatial Fourier trans- 
form of the orientationally averaged pair correlation func- 
tion. On molecular scales and microscopic times, the 
time evolution of density fluctuations of the glassforming 
Toluene and m-Toluidine is dominated by translational 
dynamics, i. e. the center-of-mass relaxational dynam- 



ics. Such a strong de Gennes narrowing effect has not 
yet been reported for the a-relaxation of low molecular- 
weight glass- forming liquids. Whether Scom(q) is the 
right input for a sophisticated mode-coupling analysis in 
the supercooled high-density state can only be decided by 
further molecular dynamics simulations pi| . Some fortu- 
itous compensation effects in the orientational fluctua- 
tions can play an important role. 

Despite the structural similarity of both molecules 
there is an important difference that might play an inter- 
esting role in the structural relaxation dynamics: the for- 
mation of hydrogen bonds in a system like m-Toluidinc 
induces clustering of a few molecules whose signature is 
a well defined prepeak around q ~ 0.5 A -1 in the static 
structure factor fn| . Unfortunately, our experiments in 
this low g-range are of too low statistical qualitiy - 
NSE experiments outside the diffraction maximum are 
extremely difficult to perfom — to decide whether this 
static heterogeneity has a counterpart in the dynamics 
or not. 
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Figure captions 



Figure 1: 

Density correlation functions $(<?,£) at 160 K for Toluene at three different wave vectors q. Lines are fits with the 
Kohlrausch function (|l]). Note, at go = 1.3A -1 , close to the first maximum of the center-of-mass static structure 
factor 5com(<?), the correlations decay much slower than for other wave vectors. 
Figure 2: 

Master curves $>(q,t) at several q for Toluene (a) and m-Toluidine (b) as obtained by rescaling times to t = 
tT a (To)/T a (T). The data fall onto a temperature-independent master curve. The dense points at long times steem from 
Fourier deconvoluted Mibemol TOF spectra S(q,u). For Toluene ®(q,i) at 200K for q = 0.9 A -1 and q = 1.3 A -1 is 
included (a). For m-Toluidine t) for q — 1.3 A -1 from a fragmentary experiment on Mibemol at 273 K is included 
in (b). Note the excellent agreement between NSE and TOF. 
Figure 3: 

(a) Wave vector dependence of the structural relaxation time (r q (T)) / (r go (T)) for Toluene obtained from fits with (|l|) 
to the neutron spin-echo master curves (cf. Fig. 2) (open symbols). Results from fits with the Fourier transformation 
of (Q) to the dynamic structure factor S(q,ui) from Mibemol at 200 K are included (closed symbols). The solid line 
is q~ 2 ScoM(q)- obtained by Monte-Carlo simulations at 160 K. The dotted lines is the experimentally observed total 
q~ 2 S m (q) at 161 K @,||- (b) The same for m-Toluidinc with q- 2 S C OM(q) at 260K and q~ 2 S m (q) at 248K. 
Figure 4: 

Wave vector dependences of the the prefactor f q obtained from fits of ([!]) to the master curves. Open squares: Toluene 
(NSE), open circles: m-Toluidine (NSE). The closed squares are plateau values from Mibemol TOF data obtained by 
time averaging &(q,t) at 140K in the range 3psec <t< 8psec and multiplied by 1.1 to match NSE data. Solid line: 
Scom(<7), dotted line S m (q) for Toluene as in Fig. 3. 
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